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First identical-pion correlations measured at RHIC energies by PHENIX are presented. 
Two analyses with separate detectors, systematics, and statistics provide consistent re- 



^ ' suits. The resulting HBT radii are moderately larger than those measured at lower en- 
O ■ ergies. The kt dependence of the Bertsch-Pratt HBT radii is also similar to previous 
. measures and is consistent with the conjecture of an expanding source. 



1. INTRODUCTION 



O 

V ! One predicted result of a first order phase transition to a quark gluon plasma in heavy 

Zj ■ ion collisions is a long lived state due to the large latent heat and corresponding reduction 
^ ■ in pressure gradients in the created system The primary method for measuring the 
space-time extent of heavy-ion collisions is two particle interferometry, which has been 
utilized in hadronic interactions from proton induced to nuclear collisions. For an incoher- 
^ I ent source (S) of identical bosons, the normalized probability of detecting two particles 
' with relative momentum q = pi — P2 and average momentum k = {pi + p2)/2 is given by 
the correlation function (C2): 

C2(g, k) = l + \ J dxS{x, k)e'^''\ (1) 

where x is the four-position, and the integral is taken over all space-time. The direct 
utilization of this technique is hindered in principle because equation |l] is uninvertible: 
the source has 7 independent variables while the measured correlation function has 6. 
Furthermore, due to limited statistics most recent heavy ion experiments have chosen 
to plot and fit the correlation function in only three of the six dimensions versus the 
Bertsch-Pratt projection of q in the out-side-long directions: 

C2(g, k) = l + \ exp{-qTout^RTout^{k) - qTside^Rxside^ik) - qLong^Rhong^ik)) (2) 

where qiong is the momentum difference in the beam direction, qrouti in the pair momen- 
tum direction, and qTside, the corresponding orthogonal direction. A benefit of such a fit 
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Figure 1. Correlation func- 
tion for 7r~ pairs from the 
EMC analysis. The one- 
dimensional projections of 
the three-dimensional cor- 
relation function is aver- 
aged over the lowest 80 
MeV/c in the other mo- 
mentum differences. 

is that, for a well-behaved, gaussian source with no position-momentum correlations or 
resonances, the lifetime is simply determined by the difference between Rrside and Rtoui- 
P'^t'^ = Rtoui' — Rrside^ ID- In such studies the k-dependence of the HBT radii is usually 
explored by repeating the fitting procedure with a subsample of the entire dataset of pairs. 

2. PHENIX DETECTOR AND ANALYSIS 

The PHENIX experiment is described in detail elsewhere . For the analysis described 
in this article we use the data collected in the summer of 2000. After all offline analysis 
cuts, the data sample was approximately 1.5 million events. 

This analysis used a subset of all of the detectors in the experiment. For vertex in- 
formation and centrality definitions we rely on a pair of Cerenkov beam counters and 
zero-degree calorimeters; the beam counters also act as the start of the time of flight 
measurements. The two drift chambers measure the deflection of particles through the 
magnetic field and, hence, their momentum, while the velocity is measured by either the 
time-of-fiight wall (TOF) or the electromagnetic calorimeter (EMC). The TOF timing 
resolution is approximately 115 ps while the EMC resolution is 700 ps. As a re- 
sult, the TOF analysis provides tt-K separation to 1.5 GeV/c; the EMC analysis, to 0.7 
GeV/c. Due to geometric and acceptance effects the TOF analysis consists of ~ 500, 000 
(~ 300, 000) pairs of identified vr"*" (vr") while the EMC analysis consists of approximately 
five times as many pairs. 

The particle identification algorithm is similar in both analyses. A pion is defined 
as being within 2a of the pion mass-squared peak but 3a away from the kaon peak. 
Backgrounds from long lived resonances, e.g. A's, are reduced by intradetector association 
cuts. Energy deposition cuts in the TOF slat or EMC tower further reduce backgrounds 
due to accidental misassociations. Systematic studies have shown that the background 
contamination to the EMC analysis is at least double that in the TOF analysis leading 
to artificially lower A values in the former. Ongoing studies of the backgrounds in the 
analyses aim to resolve these differences. However, the introduction of background from 
accidental hits or electron contamination tend to only affect the resulting A in the fits and 
do not change the measured radii. 

A number of systematic studies have been performed to ensure that the resulting cor- 
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Data Set 


Rxout (fm) 


Rrside (fm) 


Rlong (fm) 


A 


EMC 7r+7r+ 


4.4 ±0.2 


5.1 ±0.6 


5.9 ±0.4 


0.27 ± .02 


TOP 7r+7r+ 


6.2 ±0.5 


7.9 ± 1.1 


4.0 ± 1.2 


0.49 ± .05 


EMC TT-TT- 


5.1 ±0.2 


5.0 ±0.6 


5.9 ±0.4 


0.30 ± .02 


TOF TT-vr- 


5.5 ±0.5 


5.8 ± 1.5 


6.7 ±0.9 


0.49 ± .06 



Table 1 

Results of the Bertsch-Pratt fits to the identical pion pairs in the EMC and TOP analyses. 
Errors shown represent statistical uncertainties only; current systematic uncertainties are 
< 1 fm. 



relation function does not include any artificial, experimentally induced correlations from 
inefficiencies in either the detector or tracking algoithms. To remove these inefficiencies, 
pairs of particles within 2 cm of each other in the drift chamber are removed. Purther, 
pairs that share the same TOP slat or EMC cluster are also removed in both the signal 
and mixed background. Por the Coulomb correction, the source is parameterized as a 
gaussian in Rinv= VSr"^ — 6t^, as determined by an iterative procedure [Q. 

3. RESULTS 

Pigurej^ shows the Coulomb corrected correlation function for tt" pairs measured in the 
EMC analysis overlayed with the resulting fit. While the fit is performed in the full three 
dimensional space, we plot here projections of the correlation function into each of the 
standard momentum difference variables {grout, iTside, and qLong)- The results of the fit 
are shown in Table along with the results from the tt"'" analysis from the EMC and the 
7r+/7r^ analyses utilizing the TOP. The mean transverse momentum {{kf)) of the pairs is 
350 MeV/c in the TOP analysis and 340 MeV/c in the EMC analysis while the rapidity 
coverage in both analyses is centered about mid-rapidity, \y\ < .35. The mean centrality 
of all pairs in the analysis is 15% and is strongly biased towards central collisions. 

Within the current statistical and systematic error bars, the EMC and TOP analyses 
for both pion sets are consistent with one another. The results of the fits are moderately 
larger than identical measurements for Au±Au and Pb±Pb collisions at lower energies 
and comparable {kt) 0. The results do not indicate an especially large source compared 
to measurements at lower energies and the resulting naively calculated lifetime deduced 
from these measurements is consistent with zero. However, a stronger collective flow will 
shorten the effective source and lifetime measured by HBT. 

Plotted in Pig. ^ is the (kf) dependence of the radii in the EMC 7r+ (squares) and 
7r~ (triangles) analyses. The sample is split into three approximately equal subsets of 
pairs and the correlation analysis is performed on each subset. The kt bins are kt < 250 
MeV/c, 250 < h < 350 MeV/c and h > 350 MeV/c corresponding to {h) = 188 MeV/c, 
{kt) = 298 MeV/c and {kt) = 436 MeV/c respectively. Data points at identical {kt) are 
offset by 10 MeV/c for clarification. 

The radii follow trends familiar from lower energy measurements at lower beam energies 
0: Rrside has a mild, if any, dependence on the transverse momentum of the pair, while 
Rxout varies strongly with {kt). Such dependencies have been described by collective 
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Figure 2. kt dependence 
of the Bertsch-Pratt radii 
for TT^ (triangles) and tt"*" 
(squares) pairs in the EMC 
analysis. Error bars corre- 
spond to statistical errors 
only. 



motions, though temperature gradients and decays of short-lived resonances could also 
contribute to such dependencies p. 

4. CONCLUSION 

We have shown first measurements of identical pion correlations measured by PHENIX 
at RHIC. The resulting Bertsch-Pratt radii are moderately larger than those measured 
at lower energies. The kt dependence of the radii is consistent with the conjecture of a 
system with a large degree of collective motion. 

An extension of these measurements over an wider range in kt with a larger data sample 
will be available following the upcoming 2001 data collection. During this second run the 
additional PHENIX acceptance and integrated luminosity should lead to approximately a 
factor of a hundred in pair statistics. The resulting data sample will lead to much stricter 
constraints on models of the collision and allow for comparisons to identical kaon and 
proton HBT analyses. 

This work was performed under the auspices of the U.S. Department of Energy by the 
University of California, Lawrence Livermore National Laboratory under Contract No. 
W-7405-Eng-48. 
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